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ABSTRACT 

Individual pulses from pulsars have intensity-phase profiles that differ widely from 
pulse to pulse, from the average profile, and from phase to phase within a pulse. 
Widely accepted explanations do not exist for this variability or for the mechanism 
producing the radiation. The variability corresponds to the field statistics, particu- 
larly the distribution of wave field amplitudes, which are predicted by theories for 
wave growth in inhomogeneous media. This paper shows that the field statistics of the 
Vela pulsar (PSR B0833-45) are well-defined and vary as a function of pulse phase, 
evolving from Gaussian intensity statistics off-pulse to approximately power-law and 
then lognormal distributions near the pulse peak to approximately power-law and 
eventually Gaussian statistics off-pulse again. Detailed single-component fits confirm 
that the variability corresponds to lognormal statistics near the peak of the pulse pro- 
file and Gaussian intensity statistics off-pulse. The lognormal field statistics observed 
are consistent with the prediction of stochastic growth theory (SGT) for a purely 
linear system close to marginal stability. The simplest interpretations are that the 
pulsar's variability is a direct manifestation of an SGT state and the emission mech- 
anism is linear (either direct or indirect), with no evidence for nonlinear mechanisms 
like modulational instability and wave collapse which produce power-law field statis- 
tics. Stringent constraints are placed on nonlinear mechanisms: they must produce 
lognormal statistics when suitably ensemble-averaged. Field statistics are thus a pow- 
erful, potentially widely applicable tool for understanding variability and constraining 
mechanisms and source characteristics of coherent astrophysical and space emissions. 

Key words: pulsars: general; pulsars: Vela; radiation mechanisms: non-thermal; 
methods: statistical; waves; instabilities. 



1 INTRODUCTION 

Pulsars are believed to be highly magnetized neutron stars 
whose rotation causes highly nonthermal beams of radiation 
to be swept across the Earth. Most likely the radiation is pro- 
duced over the magnetic polar caps of the star, which are 
offset from the rotational poles. The radiation's high bright- 
ness temperatures require coherent emission processes such 
as plasma microinstabilities or nonlinear processes; however, 
despite many years of research, no agreement exists on which 
mechanisms dominate (Asseo 1996; Hankins 1996; Melrose 
1996; Melrose & Gedalin 1999). Proposed linear instabilities 
include: (i) Linear acceleration and maser curvature emis- 
sion (Luo & Melrose 1995; Melrose 1996), in which elec- 
trons radiate coherently while accelerating in an oscillating 
large-scale field or on curved magnetic field lines, respec- 
tively, (ii) Relativistic plasma emission (Melrose 1996; As- 
seo 1996), in which a streaming instability either directly 
generates escaping radiation near harmonics of the electron 



plasma frequency /pe or else drives localized, non-escaping 
waves near fpe that are converted into escaping harmonic 
radiation by linear mode conversion or nonlinear processes, 
(iii) A streaming instability into a new, directly escaping 
mode (Gedalin, Gruman & Melrose 2002). Possible non- 
linear mechanisms involve solitons, modulational instabili- 
ties, and strong turbulence wave collapse of intense localized 
wavepackets of waves driven near /pe (Pelletier, Sol, & Asseo 
1988, Asseo, Pelletier, & Sol 1990, Asseo 1996, Weatherall 
1997, 1998). Another possibility is an antenna mechanism, 
in which largescale low frequency wavepackets act as an- 
tennas for conversion of other higher frequency waves into 
radiation, due to acceleration of electrons in the combined 
wave fields (Pottelette, Treumann & Dubouloz 1999, Cairns 
& Robinson 2001). Standard analyses of linear and nonlinear 
growth rates suggest that numerous mechanisms are viable, 
in part due to uncertainties in the source plasma character- 
istics and location of emitting regions (e.g., above the polar 
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Figure 1. Average flux profile as a function of phase bin for the 

Vela datasct (thick line), together with three superposed individ- 
ual pulses (other lines). See Section 3 for more details. 

cap or near the light cylinder). Accordingly new approaches 
are necessary. 

Since the discovery of pulsars it has been known that 
only suitably long time averaging leads to a stable intensity 
profile. While this average profile is unique to each pulsar, 
individual pulses vary widely in intensity, often by a factor 
of 5 or more, from one phase to another in a given pulse 
and from one pulse to the next at a given phase. Illustrated 
in Figure 1 for Vela, this variability is intrinsic, clearly re- 
lated to the emission mechanism and/or propagation effects, 
and has no accepted interpretation. This paper considers the 
variability in terms of its statistics, following Cairns, John- 
ston & Das (2001), hereafter called Paper I. In general, the 
variability includes phenomena known as drifting subpulses, 
microstructures, giant pulses, and giant micropulses. Sub- 
pulses are features that drift in time across the pulse win- 
dow (Drake & Craft 1968; Manchester & Taylor 1977), while 
microstructures are concentrated features superposed on a 
subpulse that are sometimes quasiperiodic (Craft, Comella 
& Drake 1968, Kramer et al. 2002). Giant pulses (Cognard 
et al. 1996; Hankins 1996) and giant micropulses (Johnston 
et al. 2001) are very rare pulses with pulse-integrated fluxes 
~ 10 times the average. 

Analyses of field statistics, such as the distributions of 
electric field strengths or intensities, are not yet standard 
in analyses of astrophysical radiation, in contrast to Fourier 
and correlation analyses related to propagation efTects (Rick- 
ett 1977). Probable reasons include (i) the strong theoret- 
ical motivations and benefits of considering field statistics 
were not clear before the advent of stochastic growth the- 
ory (SGT) (Robinson 1992, 1995, Robinson, Cairns & Gur- 
nett 1993, Cairns, Robinson & Anderson 2000, Robinson 
& Cairns 2001, Cairns & Menietti 2001) and other theo- 
ries such as self-organized criticality [SOC] (Bak, Tang & 
Weisenfeld 1987, Bak 1996), (ii) these motivations and ben- 
efits are not widely known in the astrophysical community, 
and (iii) high time resolution, coherently dedispersed data 
for high intensity sources, whose statistics are not strongly 
contaminated by instrumental/background noise, have only 
recently become available (Johnston et al. 2001; van Straten 
et al. 2001; Kramer et al. 2002). Recent analyses of seven so- 



lar system wave phenomena show that all have well-defined 
field statistics that agree very well with the functional forms 
predicted by SGT, resolving longstanding theoretical prob- 
lems pertaining to the burstyness, widely varying fields, and 
persistence of the waves (and driving particles) for unex- 
pectedly large distances in the source. Giant pulses of some 
pulsars have power-law flux distributions (Cognard et al. 
1996), as do giant micropulses (Johnston et al. 2001), poten- 
tially being interpretable in terms of SOC (Young & Kenny 
1996), nonlinear modnlatioual processes (Weatherall 1998; 
Paper I), or driven thermal waves (Cairns et al. 2000). 



Pulsar variability and emission mechanisms can be ad- 
dressed directly by analyzing the radiation statistics. Re- 
cently we presented a preliminary analysis of field statistics 
for the Vela pulsar (Paper I), demonstrating that the ob- 
served variability is intrinsic and corresponds to lognormal 
statistics at some pulsar phases, interpreting these statistics 
and the observed variability in terms of a pure SGT system, 
so that the associated pulsar emission mechanism involves 
only linear processes. The overall goal of the present pa- 
per and its companion (Cairns et al., 2002, hereafter called 
Paper III) is to present a detailed analysis of the Vela pul- 
sar's field statistics. Specific goals of this paper are; (i) to 
survey Vela's field statistics as a function of phase, show- 
ing evolution from approximately Gaussian to power-law to 
lognormal to power-law to Gaussian distributions as the 
pulse phase varies from off-pulse to on-pulse to off-pulse 
again, (ii) present single-component Gaussian and lognor- 
mal fits to the field statistics for appropriate pulsar phases, 
(iii) demonstrate in detail that the variability corresponds 
to lognormal statistics near the peak of the pulse profile but 
Gaussian noise off-pulse, (iv) interpret the lognormal statis- 
tics in terms of SGT and linear emission mechanisms, and 
(v) discuss the results, including the constraints on nonlin- 
ear emission mechanisms. Paper III presents and interprets 
two-component fits to the observed field statistics over the 
remaining phase ranges; together these papers show that 
the Vela data are consistent with SGT applying, and with 
the emission mechanisms being purely linear, whenever the 
pulsar is detectable. 



This paper proceeds by presenting required background 
information on theories for wave statistics (Section 2) and 
the Vela dataset (Section 3). Subsequently, the observed 
field statistics are surveyed, their evolution demonstrated, 
and initial analyses and associated interpretations described 
(Section 4). Detailed single-component fits to the observed 
field distributions for specific ranges of phase are then pre- 
sented, first Gaussian intensity distributions off-pulse (Sec- 
tion 5) and then lognormal distributions near the center of 
the average pulse profile (Section 6). Brief descriptions of 
attempts to fit other distributions, specifically distribu- 
tions and Gaussians in the field, to the data for these phase 
ranges are summarized in Section 7. Theoretical interpreta- 
tions of these results are described and placed in context, 
first specifically for the Vela pulsar (Section 8) and then 
for other pulsars and astrophysical and solar system sources 
(Section 9). The conclusions are given in Section 10. 
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2 THEORIES FOR WAVE STATISTICS 

Wave growth in inhomogeneous plasmas naturally leads to 
bursty waves with widely varying fields. A qualitative ra- 
tionale for this is as follows: (i) the plasma's inhomogeneity 
leads to localized regions being favored for wave growth (e.g., 
the plasma instability depends upon plasma density, temper- 
ature etc.), (ii) the distribution function /(v) for the driving 
particles relaxes more in the growth sites (leading to larger 
wave fields there), (iii) this injects spatially and temporally 
varying fluctuations into /(v), (iv) these fluctuations cause 
burstyncss and variable waves along the particle path, and 
(v) the system evolves to a statistically steady state where 
possible. These wave-particle interactions are expected to 
drive the system towards marginal stability, where wave 
emission and damping (and related energy inflows and out- 
flows) are balanced, time- and volume-averaged. One likely 
condition for reaching a steady state is that the unstable 
features in /(v) be able to reform, at least partially; e.g., 
due to fast particles outrunning slow particles in a beam 
system (Robinson et al. 1993). Different theories for wave 
growth are then characterized by different degrees of inter- 
action between the waves, driving particles, and background 
plasma, leading to different wave statistics. Since these the- 
ories are not well known in the astrophysical literature, their 
salient features are summarized here in some detail. 

SOC (Bak et al. 1987; Bak 1996) is relevant to systems 
with fully self-consistent interactions between the waves, 
driving particles and background plasmas and with no pre- 
ferred distance or time scales. It predicts power-law statistics 
for the distribution of energy releases (e.g., wave energies), 
so that the distribution P{I) of wave intensities I obeys 
(where I for the wave electric field E) 

P(7)(x7-". (1) 

Typically the indices a approximately equal 1, with a range 
w 0.5 — 2. Proposed examples include tokamak turbulence 

(Carreras et al. 1996) and solar flares (Lu & Hamilton 1991). 
Also, Jovian "S" bursts have a power-law distribution for the 
radiation flux (Queinnec & Zarka 2001) and so may be an 
SOC system. 

SGT (Robinson 1992, 19951 Robinson et al., 1993, 

Cairns et al. 2000, Robinson & Cairns 2001, Cairns & Mc- 
nietti 2001) treats systems in which self-consistent wave- 
particle interactions occur in an independent, spatially in- 
homogeneous medium. The medium and wave-particle inter- 
actions then determine the relevant distance and time scales. 
In SGT the system evolves to a state in which (i) /(v) is 
close to time- and volume- averaged marginal stability but 
with fluctuations about this state which cause (ii) the wave 
gain G to be stochastic variable. Relevant definitions are 

£;'(t) = Boe'^**' , (2) 
with G = 2ln{E{t)/Eo) for a reference field Eq, and 

G{t)= I dt^it) (3) 

«/ — oo 

where 7(t) is the (amplitude) growth rate of the waves. SGT 
is then a natural theory for bursty waves with widely varying 
fields (due to log E oc G being a Gaussian random variable) 
that persist with the driving particles for unexpectedly large 
distances and/or times (due to the closeness to marginal 
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Figure 2. Predictions for the P(log E) distribution for pure SGT 
(solid), SGT with a nonlineaj: decay process (dashed line), and 
SGT with wave collapse (dotted line). 

stability). Moreover, SGT should be widely applicable since 
there is a natural qualitative route to an SGT state: writing 
the time integral in (3) as 

/ dt 7(t) = EiAGi , (4) 

«/ — oo 

then provided that sufficient fluctuations AGi [associated 
with fluctuations in /(v)] pass though a growth site during 
the characteristic time waves grow there, the Central Limit 
Theorem (CLT) predicts that G will have Gaussian statistics 
irrespective of the detailed statistics of AG,. This means 
that P{logE) cx P{G) obeys 

P(log£;) = {V2^a)-' ^ 

where fj, = (log E) and a are the average and standard de- 
viation of logE, respectively, and log = log^Q. That is, pure 
SGT predicts lognormal statistics for the field (Figure 2). 
In recent years pure SGT has been shown to be widely ap- 
plicable (Robinson et al. 1993, Cairns and Robinson 1997, 
1999, Cairns et al. 2000, Cairns & Grubits 2001, Cairns & 
Menietti 2001, Paper I). 

SGT can coexist at moderate E with nonlinear pro- 
cesses active at high E, causing characteristic modifica- 
tions to the predicted lognormal statistics at fields close to 
and above the threshold field Ec for the active nonlinear 
process. The presence or absence of active nonlinear pro- 
cesses can thus be constrained using field statistics, as used 
successfully in several applications. Two types of modifica- 
tion are known, as illustrated in Figure 2. First, when a 
three- wave decay process takes energy out of the measured 
waves, the P{\ogE) distribution is cutoff at E ^ Ec with 
known analytic form (Robinson et al. 1993; Cairns & Me- 
nietti 2001; Cairns & Grubits 2001). Second, the nonlinear 
process of wave collapse (Robinson 1997, Weatherall 1997, 
1998), due to a self- focusing or modulational instability of 
a wavepacket, leads to a power-law tail at 7? ~ 75c. These 
indices are typically large: collapse of Langmuir waves in 
non-rclativistic electron-proton plasmas leads to 

P{E) cx E'^P{log E) (x (6) 

with a varying between 4 and 7, depending on the dimen- 
sionality (2-D to 3-D) and shape (isotropic versus oblate ver- 
sus prolate) of the wave packets and whether the collapse is 
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subsonic or supersonic (Robinson & Newman 1990; Robin- 
son 1997). The characteristics of collapse, or even whether 
it occurs, are not so clear for the strongly magnetized and 
relativistic electron-positron plasmas relevant to pulsar mag- 
netosphcros: while a purely electrostatic analysis (Pellctier 
et al. 1988) suggests that solitons are stable to modulational 
instability and do not collapse, electromagnetic simulations 
suggest instead that the solitons are subject to modula- 
tional instability and do collapse (Weatherall 1997, 1998). 
Although the field statistics of these collapse events are not 
known, the qualitative behaviour is very similar to that for 
collapse in electron-proton simulations (Weatherall 1997, 
1998). Accordingly, it is presumed hereafter that, if mod- 
ulational instability occurs, then it leads to wave collapse 
and power-law field statistics described by (6) with indices 
similar to those given above for electron-proton simulations. 

SGT also describes the statistics of waves driven by a 
linear instability from the thermal level Eth (Robinson 1995; 
Cairns et al. 2000), being essentially a power- law with low 
index at fields ^ Eth, while thermal waves have a field dis- 
tribution given by the product of a power-law with a Gaus- 
sian in E^. In contrast, elementary burst theory [EBT] pre- 
dicts exponential statistics (Melrose & Dulk 1982, Robin- 
son, Smith & Winglee 1996) while the formerly standard 
theory for wave growth (Krall & Trivelpiece 1973; Melrose 
1986), that of homogeneous exponential growth with con- 
stant growth rate until saturated at high fields by a nonlinear 
process or quasilinear relaxation (uniform secular growth), 
predicts P(log E) should be uniform (fiat) below the nonlin- 
ear level (Robinson et al. 1993; Cairns & Robinson 1999). 
While EBT appears consistent with solar microwave spike 
bursts (Robinson et al. 1996; Isliker & Benz 2001), the statis- 
tics of all systems published to date are inconsistent with 
uniform secular growth. 

The final model considered here for field statistics and 
variability is that of Gaussian statistics in I, i.e., 

P(7) = (V2^cT/)-' g-(/-/o)V2<^? ^ (7) 

where Jo — (!) and aj arc the average and standard devia- 
tion of /. This is the standard expectation for background 
noise, for instance due to superposition of multiple random 
signals, so that the sky background, measurement noise, 
and sources with multiple unresolved sources should have 
P{I) distributions with Gaussian statistics. Scattering of 
waves between the source and observer, due to refraction 
by density irregularities, can also lead to Gaussian P{I), as 
sketched next. The Central Limit Theorem for a monochro- 
matic (real, scalar transverse) field E predicts Gaussian field 
statistics if sufficiently many ray paths contribute, whence 
the distribution P{I) should then have an exponential dis- 
tribution (RatclifFe 1956; Rickett 1977), corresponding to a 
distribution with 2 degrees of freedom. This distribution 
is well approximated by (7) if the scattering is weak and/or 
(Ratcliffe 1956). Intuitively, inclusion of finite band- 
width effects, due to multiple independent monochromatic 
signals contributing in a receiver bandwidth, should also 
lead towards Gaussian intensity statistics. 

Relationships between the P{I), P{E) and P(logi5) 
distributions follow from their normalization conditions, 
e.g., J dX P{X) = 1, and differentials. Accordingly 

P(log E) = In 10 S P(£) = 2 In 10 E^ P{I) , (8) 



allowing (1), (5), and (6) to be related quickly. In summary, 
the field statistics predicted by the above theories are known 
and can be compared robustly with observational data to 
(i) determine the relevance of particular theories for wave 
growth, thereby constraining the physics of the wave growth 
and source region and the causes for the source's variabil- 
ity, (ii) identify the presence or absence of nonlinear pro- 
cesses actively participating in the emission processes, and 
(iii) constrain the importance of scattering. 



3 DATASET 

The dataset consists of 20, 085 contiguous intensity pulses of 
the Vela pulsar, spanning approximately 30 minutes, mea- 
sured at 1413 MHz by the Parkes radio telescope and de- 
scribed in detail elsewhere (Johnston et al. 2001). Only 
a brief summary is given hero. The measured quantities 
are the voltages (proportional to the received electric field 
strength) on two orthogonal feeds, which are frequency- 
downconverted, the DC offset removed with a time constant 
tdc long compared with the pulsar period (« 88 ms), and 
fed into a modified version of the Caltech and Princeton 
backend systems (Jenet et al. 1997; Stairs et al. 2000). The 
backend system then samples the datastream in quadrature 
with 2-bit accuracy at 20 MHz and writes it to digital linear 
tape. Offiine, dispersion effects are removed using coherent 
dedispersion techniques, the four Stokes parameters are re- 
covered, instrumental polarization effects removed, and a 
flux calibrator (Hydra A) used to write the output data in 
terms of flux densities F averaged over the backend's 20 
MHz bandwidth. The resulting dataset has 2048 phase bins 
per pulse period, each of ~ 44 fis length, which is compara- 
ble to the scatter-broadening time at this frequency. About 
1800 2-bit samples are thus used to calculate the flux in 
each phase bin every pulse. Arguments against this digitiz- 
ing procedure significantly modifying the true field statistics 
are presented in section 9. 

The resulting brightness temperature Tbr{4>) at a given 
pulse phase </>, proportional to F and ideally to the antenna's 
electric fleld squared, may be written as 

Tbr{(t>) = Tree + Tgnr + Tfcacfc + Tpsr{<i>) , (9) 

where Tree is the receiver (noise) temperature, Tsnr is the 
temperature of Vela's supernova remnant, Tbacfe is the back- 
ground temperature due to galactic synchrotron radiation, 
the Sun, and other sources, and Tpsr is the pulsar's bright- 
ness temperature. Since the antenna points at the same lo- 
cation of the sky for the observing interval, removing the 
DC offset should eliminate Tsnr and Tback from the final 
dataset (provided they do not vary on times ^ tdc)- Sim- 
ilarly, slow drifts in Tree will be removed, although rapid 
fluctuations corresponding to the receiver's thermal noise 
persist and naturally dominate the observed field statistics 
at off-pulse phases (see Section 5 below). 

While removal of the slowly- varying DC offset is an ad- 
vantage in terms of observing weak sources above receiver 
and sky backgrounds, it is also a disadvantage for analy- 
ses of field statistics. There arc at least three reasons. First, 
since E^ (x I oc Tpsr, neglecting the time-steady, phase- 
averaged part of Tpsr {{Tpsr}) affects the field scale non- 
linearly for fields E such that Tbr{4>) is within a factor of 
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a few times (Tpsr), thereby potentially modifying the field 
statistics there. Second, the logarithim function, which fea- 
tures prominently in SGT and the observed field statistics 
shown below, is undefined for negative fluxes (which occur 
in the dataset due to removal of the DC offset and duo to 
thermal fluctuations in the receiver noise). Third, measure- 
ment of non-zero values of (Tpsr) would permit more de- 
tailed investigation of the source physics, particularly if mul- 
tiple wave populations are superposed, and perhaps allow 
competing physical models to be tested and distinguished. 
Possible origins for physically-significant, non-zero (Tpsr) in- 
clude synchrotron emission from the pulsar magnetosphere 
and coherently-produced radiation that undergoes scatter- 
ing and diffusion as it propagates from its source, perhaps 
even changing the phase at which it is observed. Since the 
values (Tpsr-) are not available for these data, these issues 
are addressed by adding a constant offset flux 

7;yy = 1250 mJy (10) 

to each sample in the dataset. Clearly, when fltting Gaussian 

intensity distributions (7) to the data, only relative varia/- 
tions in the centroid from 1250 mjy arc significant. 

Comparisons between the data and Section 2's theories 
are most easily accomplished by recasting the data in terms 
of fields and intensities. This is not inappropriate anyway, 
since the radiation fields incident on the antenna were first 
detected as voltages and the final calibrated fluxes arc aver- 
aged over the detector's 20 MHz bandwidth. The calibrated 
intensity / then equals the flux multiplied by the bandwidth 
and the field E is proportional to 7^^^. Accordingly, rather 
than converting into SI units, the analyses below are cast in 
terms of variables E' and /', defined by 

E' = {F + 1250mJy)^''2 oc E , (11) 
I' = F + 1250mJy oc I , (12) 

which are related directly to the calibrated fluxes F and 
whose units arc in (mjy)^'''^ and mJy, respectively. Accord- 
ingly, E' and /' differ by multiplicative constants from E 
and I (and ideally to the incident field and intensity also), 
respectively, which correspond to constant scale factors or 
offsets along the abscissa axis in the linear and logarithmic 
plots of field statistics below. 

Figure 1 shows the average pulse profile for relevant 
phase bins in mJy (Paper I), together with three superposed 
pulses that illustrate the variability. Note that the noise level 
is very low compared with many earlier analyses, allowing 
investigation of the intrinsic field statistics. Additionally, the 
pulses are sampled rapidly in time, reaching the limits im- 
posed by scatter-broadening. This permits detailed investi- 
gation of the field statistics as a function of pulsar phase. 
The significant variations in field statistics with phase shown 
below suggest that rapid sampling is important if accurate 
analyses of field statistics is desired. 



4 IDENTIFICATION OF REGIONS 

The off-pulse phase bins considered (phases 390 — 420 and 
600 — 630, cf. Figure 1) have average intensities I' w I'^ff = 
1250 mJy and maximum values w 8000 mJy. Only the on- 
pulse bins have I' ~ 10"* Jy. Figure 3 surveys the whole 
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Figure 3. (Top) The maximum value of log E'((j>) (circles) and 
three estimates for the averages /i(^) = (log-B'}(^), calculated 
using three I' thresholds (solid - 10* mJy, dashed - 2.5 x 10^ 
mJy, and dotted - mJy), as functions of phase cp for the Vela 
dataset. (Bottom) Three estimates for (j{cj>), using the line styles 
and intensity thresholds for (Top) , as functions of phase. 



dataset as a function of (p: the top panel shows the maxi- 
mum value of logiS' for each 4>, jj,max{4') (circles), and three 
averages for ju(0) = (logS')(</)) (solid, dashed, and dotted 
lines), corresponding to averaging log E' samples above spec- 
ified intensity thresholds (/' > 10**, 2500, and mJy, respec- 
tively). The Figure's lower panel shows three estimates for 
(j{(j>), which is the standard deviation of logJS', for these 
same phases and thresholds. 

The sharp, very localized peak in /i((?!)) and (j{<j>) for a 
threshold of lO** .Jy at phases 429 — 433, and its coincidence 
with the localized peak in fimax, demonstrates the appear- 
ance of a strongly phase-localized, intense component of the 
pulsar's output. This corresponds to the giant micropulses 
identified previously (Johnston et al. 2001). Otherwise Fig- 
ure 3 shows smooth variations in /i, ^max and a with phase, 
including the "bump" identified for bins 550 — 600 by John- 
ston et al. (2001). The effects of the different thresholds are 
not important for phases « 460 — 540, where the curves 
for the three thresholds agree well, corresponding to phases 
where At(0) ~ 0.5 log \/E>ff and the pulsar fields are typi- 
cally well above the average value. However, the peaking of 
II and limax at different phases (w 480 versus « 450) points 
towards different field statistics in these regions, as demon- 
strated in detail next, and to evolution in the number and 
characteristics of active components for the source. 

Figure 4 shows the evolution of the field distribution 
PilogE') as a function of pulsar phase. At each phase this 
distribution is calculated by binning the intensity time series 
in logE' and normalizing so that J dilog E')P {log E') = 1. 

The P (logs') distributions at phases ~ 428 and ~ 600 cor- 
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respond to a Gaussian intensity distribution, being inter- 
pretable as Gaussian noise, as shown in detail in Section 
5 below. However, for <^ ~ 429 first a localized high-/ tail 
appears and disappears [the giant micropulses (Johnston et 
al. 2001)], after which an approximately power- law compo- 
nent appears at about phase 435 and extends to increas- 
ingly high E' until about phase 455. After this the entire 
distribution moves towards higher E' , broadens, and the 
high-i?' power-law slope becomes approximately parabolic 
in log P (log -E') — logiJ' space at phase 465, correspond- 
ing to a lognormal distribution. The peak of this lognormal 
P(log-B') distribution continues to move to higher E' until 
about phase 480 but becomes narrower, corresponding to 
higher /i and smaller a in (5). Subsequently, the distribu- 
tions continue to be approximately parabolic but move to 
lower fj, and a. Eventually the effects of Gaussian intensity 
noise become evident, giving rise to a non-parabolic tail at 
low E' for phases ~ 540. In the phase range « 545 — 600 
a power-law tail at high E' develops and then disappears, 
corresponding approximately to the "bump" (Johnston et 
al. 2001) in /x in Figure 3, after which the field distribution 
recovers the form in the off-pulse bins below phase 428. 

Figure 4 thus directly demonstrates that the Vela pulsar 
has different field statistics in neighboring ranges of phases; 
put another way, the field statistics vary with across the 
source. This immediately shows that field statistics can be 
used to probe the source physics. 

The qualitative identifications given above are elabo- 
rated next by considering the quantities 



(13) 
(14) 



the first of which was considered previously, for example, by 
Johnston et al. (2001). A Gaussian intensity distribution at 
phase (f) is expected to have Rjf (</>) ~ 5, since almost all its 
samples should be within ~ 5(7// (i/i) of the mean A 
similar statement follows for a lognormal field distribution, 
which should have Rfi{4') ~ 5. Figure 5 shows these quan- 
tities at different phases for the three intensity thresholds 
used in Figure 3. The Rji{(f)) values in Figure 5's top panel 
show that Vela's field statistics are close to Gaussian in 7' 
only for phases ~ 428 and ~ 610. The ii^ values in Figure 
5's bottom panel show that the field statistics are close to 
lognormal at phases ~ 460 and ~ 540, with particularly close 
agreement evident in the range 460 — 510 where the traces 
for all three threshold are very close together. 

Figures 3 to 5 are thus all qualitatively consistent with 
the following identifications: phases ~ 428 and ~ 600 have 
Gaussian intensity statistics (Paper I), phases 429 — 433 
have non-Gaussian and non-lognormal statistics that cor- 
respond to giant micropulses (Johnston et al. 2001; Kramer 
et al. 2002), phases with 435 - 455 and 545 - 600 have an 
approximately power-law character at largo /' and Iog_B'; 
phases 460 — 540 have lognormal field statistics (Papers I and 
III). These identifications are quantified next, using single- 
component fits, for the Gaussian-7' and lognormal regimes, 
refining and extending the analysis of Paper I. The approx- 
imately power-law domains are termed the "transition re- 
gion" in Paper III, since these distributions are best inter- 
preted in terms of vector convolution of a Gaussian and a 
lognormal component and not in terms of intrinsic power- 
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Figure 4. Distributions P(log£"') for chosen piiascs 4>, ranging 
from off-pulse (~ 430) to pulse-center (e.g., 470 and 490) to off- 
pulse again (~ 600). 



law statistics. Two-component fits for the transition region 
and lognormal domains are detailed in Paper III. 



5 GAUSSIAN NOISE REGION 

The field statistics at phases ~ 428 and ~ 600 are now shown 
in detail to be Gaussian in the intensity /', as expected for 
instrumental noise, background "noise" formed by incoher- 
ent superposition of a large number of small signals, and/or 
scattering. Figure 6 shows the P(/') distribution observed 
for phases 391 — 400, calculated by binning the data into 
linear intensity bins and normalizing, together with the cor- 
responding fit to (7), obtained using the Amoeba algorithm 
to minimize (Press et al. 1986). Restricting the fit to 
intensity bins with at least 100 samples, the fit parame- 
ters are (/') = 1220 mJy (agreeing with I'^ff to within less 
than the 100 mJy bin width), aj, = 1420 mJy, = 67 for 
Ndof = 47 degrees of freedom, and a significance probabil- 
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Figure 6. Comparison of the observed intensity distribution 

P(I') for phases 391 — 399 (circle symbols with error bars given 
by counting statistics) with the best fit to the prediction (7) for 
Gaussian statistics (solid line). Fits are performed where the dot- 
ted line shows that bins have > 100 counts and to the right of 
the vertical dashed line at /' = 200 mJy. 



ity -P(X^) = 0.03. Clearly the fit (7) agrees well with the 
observed data and has reasonable statistical significance. 

At high I' ~ 7000 mJy the observed distribution hints 
at deviations from the Gaussian form (7). The simplest ex- 
planation is that these high- 7' samples correspond to pulsar 
emissions that slightly exceed the background. Put another 
way, the phase range 391-400 is not entirely ofi'-pulse. 

Lognormal fits to Figure 6's data are clearly inferior 
(not shown). This is expected based on the strong depen- 
dence of the R^{<t>) curves in Figure 5 on the intensity 
thresholds used, in contrast to the lack of variation of the 
Rji {(j)) curves on these thresholds. Thus these data are con- 
sistent with Gaussian intensity distributions and not log- 
normal field (or intensity) distributions. Analyses of other 
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Figure 7. Comparison of the intensity distribution P{I') ob- 
served at phase 490 with the best fit (7) for Gaussian statistics 
(solid line). The circle symbols, error bars, dotted line, and ver- 
tical dashed line axe as in Figure 6. 



phases in the domain stated above yield similar fit parame- 
ters and statistics. This is expected on the basis of the very 
similar P(logS') distributions shown in Figure 4. 



6 PURE LOGNORMAL REGION 

Single-component lognormal fits are presented in section 6.1 

for the on-pulse phase domain 470 — 540. Suggestions that 
a second component may also contribute are discussed in 
section 6.2. Two-component fits, which lead to very good 

agreement between observation and theory over the entire 
phase range, are presented in Paper III. 



6.1 Single component fits 

Figure 7 shows the P{I') distribution and its best Gaussian 

fit for phase 490, close to but after the peak in the pulsar's 
average profile (Figure 1). The fit clearly fails at both low 
and high /', entirely missing the long tail at large /', and has 
= 301 for Ndof = 53 and P(x^) < lO"^*^ for the fitted 
range of I' (dotted horizontal line). The variability at this 
phase is thus not described by Gaussian intensity statistics. 

In contrast. Figure 8 shows that the P(logi?') distribu- 
tion for phase 490 is well fitted by the pure SGT prediction 
(5): for bins with > 100 counts (dotted line) and E' > 10^ 
(mJy)^/2 (/' > 10'' mJy, which is Qaj above {/') in Figure 
6), the fit parameters are /i = 2.3, a = 0.096, = 27 for 
Naof = 19 and P(x^) = 0.12. Similarly the Kolmogorov- 
Smirnov test (Press et al. 1986) yields a significance prob- 
ability of 0.47. This fit is strongly statistically significant: 
pulsar variability at this phase is lognormally distributed 
and quantitatively consistent with the theoretical form (5) 
predicted for pure SGT. Note that the fit matches the data 
well even outside the fitted range of fields (vertical dashed 
line and horizontal dotted line), although the effects of the 
noise background become increasingly evident at fields < 80 
(mjy)^/^. It is worth emphasizing that the P{I') distribu- 
tion in Figure 7 for 3 x 10* ^ /' S 1.7 x 10^ mJy does 
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Figure 8. Comparison of the distribution P(log E) observed at 
phase 490 witli the best fit (5) for lognormal statistics (soUd line). 
The format is identical to Figure 6 except that the vertical dashed 
Une is now at E' = 100 mjyi/^. 
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Figure 10. Comparison of the distribution 1'{X) for all data 
in phases 485-495 (plus signs), calculated using the procedure 
described in the text, with the pure SGT prediction (solid line). 
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Figure 9. Comparison of the distribution P(logi?') observed at 
phase 510, in the same format as Figure 8, with the best fit to 
the prediction (5) for lognormal statistics (solid line). 



not have an intrinsically power-law form but is instead best 
represented in terms of a lognormal form (Figure 8). 

Figure 9 presents the observed P(log E') distribution 
and fit to the pure SGT prediction (5) for phase 510 . Again 
the observed distribution is well fitted by pure SGT with 
good statistical significance. For fields E' > 130 (mjy)^/^ 
the fit parameters are fj, = 2.1±0.1, a — O.lOiO.Ol, = 13 
for Ndof = 6 with -P(x^) = 0.04. The statistical significance 
changes somewhat with the fitting threshold in E': requir- 
ing E' > 10^ (mJy)i/2 yields /x = 2.1 ± 0.1, a = 0.11 ± 0.1, 
= 42, Ndof = 12 and P(x^) = 3 x 10"^, while requiring 
E' > 150 (mJy)^/^ yields ^ = 2.1 ± 0.1 and a = 0.10 ± 0.01 
but x^ = 7 for Ndof = 3 and P(x^) = 0.07. These varying 
statistical significauccs correspond to the varying contribu- 
tion of the noise background to the observed distribution, 
due to the relatively weak pulsar fields at this phase and to 
fi approaching the background noise level. In each case the 
form of the distribution at high log E' is very well fitted by 
the SGT prediction. 

Results similar to Figure 8-9 are found for phases 



470 — 540, although the statistical significance varies. Rather 
than showing more results for individual phases, Figure 10 
shows the distribution P{X) observed for phases 485 — 495 
(Paper I). Here X = [log E' — ij,{4>)]/ a{4>) is the field variable 
resulting from detrending variations in n and </> with phase (p. 
Comparison with (5) shows that purely linear SGT predicts 
the P{X) distribution to be Gaussian with zero mean and 
unit standard deviation (Cairns & Robinson 1999). Figure 
10 demonstrates that the observations and SGT prediction 
agree well. The statistical significance is also good, with the 
Kolmogorov-Smirnov test yielding a significance probability 
of 0.1%. 

One test of the results of Figures 8 - 10 is to compare 
the fits for fj, and a with Figure 3's values calculated directly 
from the data. Figure 11 shows that these values agree well, 

confirming that the fits arc accurate and the lognormal com- 
ponent dominates the field statistics at these phases. 



6.2 Initial evidence for a second component 

Figures 8-10 show that the pure lognormal form (5) fits the 
data very well above the peak in the P(logE') distribution 
but less well at low E' . Specifically the observed distribu- 
tion lies above the fit to (5) at low E' , suggesting a contri- 
bution there from a second group (or component) of waves. 
This is not unexpected in view of the Gaussian statistics 
found at off-pulse phases (e.g., Figure 6), which if it corre- 
sponds to measurement noise, scattered radiation, or con- 
tributions from multiple unresolved, incoherently summed 
sources, might be expected to contribute at all phases. 

In Paper III it is shown that the off-pulse Gaussian 
"noise" in Figure 6 indeed persists into the on-pulso bins 
and perhaps evolves into a second lognormal component. 
Performing two-component fits to the observed P(log E') 
distributions then leads to very good agreement between 
observation and theory over almost the entire range of E' 
for all on-pulso bins. In particular, the approximately power- 
law P{logE') distributions in Figure 4 correspond to vector 
convolution of the Gaussian component with an emerging 
lognormal component, while the clearly lognormal distribu- 
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Figure 11. Comparison of fi and a (circle symbols) obtained by 
fitting tlie P(logE') distributions to the SGT prediction (5) with 
the values obtained directly from the data set for /' > 2500 mjy 
in Figure 3. 

tions observed at phases 460 - 540 are best modelled in 
terms of convolution of two lognormal components. 



7 ATTEMPTS TO FIT OTHER 
DISTRIBUTION FUNCTIONS 

Attempts were made to assess the uniqueness of the Gaus- 
sian and lognormal fits presented heretofore, by also consid- 
ering several other fitting functions: a Gaussian in the field 

(rather than the intensity), the SGT prediction for a thermal 
wave distribution (Robinson 1995; Cairns et al. 2000), and 
a distribution in /'. Both the first and last of these are 
related to scattering theory (Ratcliffe 1956; Rickett 1977), 
the first being predicted for a scattered monochromatic real 
field and the latter (with 2 degrees of freedom) for the in- 
tensity of a monochromatic scattered field. In more detail, 
a distribution with n degrees of freedom corresponds to 
the distribution expected from summing a limited number 
n of Gaussian-distributed variables and it is defined by 

Px(n, /') = {l'/h)''^''-'e-''^'"^ /2"/''r(n/2) . (15) 

Here Ji is related to the average (/') by Ji = {I')/{n -\- 2). 
In the limit n ^ oo a, distribution tends to a Gaussian. 

Applying these three alternatives led to unsatisfactory 
results for all off-pulse phase bins (not shown), both in ab- 
solute terms and relative to the Gaussian intensity fits do- 
scribed in Section 5. This implies that the off-pulse data are 
best modelled in terms of Gaussian intensity statistics. 

The results of these fits for on-pulse bins in the range 
470 — 540 were also unsatisfactory, both in absolute terms 
and relative to the lognormal fits in Section 6. Specifically: 
(i) These data were not well fitted as being Gaussian dis- 
tributed in S' or as a thermal SGT distribution for any 
phase, (ii) Fits to the chi-squared distribution (15) invari- 




I (mJy) 

Figure 12. Comparison of the distribution P{I') observed at 

phase 490 witli the best fit to prediction (15) for a chi-squared 
distribution, in the same format as Figure 6. 



ably failed at high \ogE' where the lognormal fits were 
clearly superior, due to an inability to fit the high-7' tail, 
(ill) While the data could be well fitted by (15) for some 
phase bins, but not all, the fit parameters n and h var- 
ied widely rapidly from phase bin to phase bin. (iv) Both 
the statistical significance of the fits and the size of the do- 
mains in \ogE' for which good agreement existed between 
the P(log_E') data and fitting curves were generally consid- 
erably larger for the lognormal form (5) than for the chi- 
squared distribution (15). These statements are illustrated 
in Figure 12 for phase 490: for this phase (but not all) there 
is very good agreement at low /', as well as the typical inabil- 
ity to fit the high-7' tail. Moreover, the fit to (15) for phase 
490 and bins with in excess of 100 samples has n = 12.3, 
Ji = 3560, Naof = 54, = 93 and P(x^) = 8 x lO""*, cf. 
the results in section 6, while the best fit for phase 510 has 
n = 8.3, 7i = 1860 mJy, and P(x^) = 0.02. Accordingly, it 
is concluded that the on-pulsc data arc not well described in 
terms of Gaussian field statistics, thermal SGT, or a dis- 
tribution in 7', but are instead best described (for the four 
fitting functions attempted here) as lognormally distributed. 



8 THEORETICAL IMPLICATIONS FOR THE 
VELA PULSAR 

In general, the observed field statistics are determined by in- 
trinsic field statistics produced by relevant generation mech- 
anisms, possible spatial variations across the source, and 
propagation effects. Propagation effects are discounted in 
the present context, since weak or moderate scattering of 
radiation by density fluctuations is typically predicted to 
cause Gaussian or exponential statistics of 7' (Ratcliffe 1956; 
Rickett 1977), not lognormal statistics. In this dataset, then. 
Vela's variability is not due to scattering by density irreg- 
ularities. In general, though, scattering is expected to play 
some role in determining the field statistics of pulsars. The 
simplest interpretation is therefore adopted now, that the 
field statistics observed at each phase are intrinsic and are 
not significantly affected by spatial variations in the source. 
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The possibility of two generation mechanisms being active 
at a given phase is deferred to Paper III. 

For the phase range 470 - 540 considered above, where 
the Vela pulsar's average pulse profile is well above the noise, 
the field distributions do not have power-law tails or nonlin- 
ear cutoffs but instead are consistent with lognormal statis- 
tics. Accordingly, these data are consistent with the observed 
variability being a direct manifestation of a simple SGT 
state, with no evidence for SOC or uniform secular growth. 
Put another way, Vela's variability is due to lognormal field 
statistics and is consistent with the emission taking place in 
a source plasma that is in a simple SGT state. 

The simplest interpretation is, then, that the absence 
of a power-law tail or cutoff in the P{logE') distributions 
for these phases is inconsistent with nonlinear processes 
(e.g., Pelletier et al. 1988, Asseo et al. 1990, Asseo 1996, 
Weatherall 1997, 1998) such as wave collapse, modulational 
instability, and three-wave decay processes playing a sig- 
nificant role. Instead, the lognormal statistics observed are 
consistent with linear emission mechanisms, implying that 
a plasma instability in a SGT state either directly generates 
the radiation or else generates non-escaping waves that are 
transformed into escaping radiation by linear processes (e.g., 
mode conversion) alone. Moreover, these results are consis- 
tent with recent conclusions (Melrose & Gedalin 1999) that 
linear mechanisms are favoured on theoretical grounds and 
are made more topical by recent proposals of new linear 
emission mechanisms (Gedalin et al. 2002). 

At a deeper level, however, nonlinear mechanisms can- 
not be ruled out entirely. Instead, stringent conditions are 
imposed: a nonlinear mechanism is viable only if it produces 
lognormal statistics when averaged over a suitable ensemble 
of nonlinear wavepackets or structures, with no evidence for 
a power-law tail or cutoff. Since existing simulations of col- 
lapse yield power-law statistics when ensemble-averaged in 
a homogeneous system (Robinson & Newman 1990; Robin- 
son 1997), mecliaiiisiiis involving collapse arc not likely to 
be viable. 

From the definitions of jj, and a and the intensity de- 
creasing with distance R as R~^, it is easy to show that 
ij{R) is constant and that fi{R) = ^(-Ro) — log(-R/i?o), where 
/i(-Ro) is the value at the source's edge {R = Rq). Taking 
the values fi = 2.0 and ct = 0.1 to be representative of these 
phases, the distance R — 350 pc for Vela, and the value 
Ro = 30 m, yields fi{R()) ~ 20. The value i?o = 30 m results 
from assuming that the overall source is annular, with ra- 
dius equal to the neutron star radius « 10 km, and dividing 
by the 2048 phase bins used for Vela. Accordingly, the ratio 
Ho/cr « 200 in the source. The values fj,o and a will constrain 
future theoretical models for why SGT applies, similar to ex- 
isting models for solar system phenomena (Robinson et al. 
1993; Cairns & Robinson 1999; Cairns & Menietti 2001). 



9 DISCUSSION 

The foregoing analyses are the first detailed applications of 
SGT to propagating EM radiation and, simultaneously, to 
extra-solar system sources (see also Paper I). Their success 
implies that radiation statistics are an underappreciated and 
potentially very powerful tool in astrophysics (and space 
physics), and suggests that SGT may well be widely appli- 



cable to coherent astrophysical sources. Of course, SGT is 
not likely applicable to all sources or indeed to all compo- 
nents of pulsar emissions, as discussed further in Paper III. 

The importance of field statistics as a probe of source 
physics is shown above in terms of emission mechanisms. 
This is also shown empirically in terms of source structure 
by the demonstrated evolution in field statistics as a function 
of Vela's phase, in particular by the field statistics ranging 
from Gaussian in I off-pulse to approximately power-law, 
lognormal, and power-law on-pulse and to Gaussian in I off- 
pulse again. This is the first detailed investigation of field 
statistics for a pulsar and the first demonstrated evolution 
in field statistics across a source. Quantitative discussion of 
how the fit parameters fi, a etc. of the field statistics vary 
across Vela are deferred (Paper III), as is discussion of how 
representative the Vela results are to other pulsars. 

It is appropriate to discuss the possibility that the in- 
trinsic field statistics are modified by the digitization, sum- 
mation, and coherent dedispersion procedures associated 
with the telescope backend system and consequent process- 
ing of the dataset (see section 3). Four observational results 
that argue against this possibility are the following: (i) the 
off-pulse statistics are Gaussian in /', as expected for re- 
ceiver thermal noise; (ii) the on-pulse field statistics evolve 
smoothly with phase and the pulse profile (see section 4) 
and have well-defined functional forms; (iii) power-law field 
statistics are obtained from this dataset for Vela's giant mi- 
cropulses (Kramer et al. 2002) with indices of order those 
observed for giant pulses from other pulsars (Cognard et al. 
1996; Johnston & Romani 2002; Romani & Johnston 2001; 
Cairns 2002); (iv) the field distributions observed on-pulse 
can be fitted very well with lognormal functions (Paper I) 
and combinations of lognormal and Gaussian functions (Pa- 
per III), and can be interpreted theoretically in terms of 
existing theories for wave growth in plasmas; None of these 
results are expected a priori if measurement and analysis 
techniques have modified the intrinsic statistics significantly. 
Moreover, result (iii) is consistent either with the backend 
system and processing system not modifying the intrinsic 
field statistics or else with all analyses of giant pulses being 
similarly flawed. Thus, while modification of intrinsic field 
statistics by the measuring process is possible and must be 
kept in mind, the above points are strong arguments against 
Vela's intrinsic field statistics being significantly modified for 
this dataset. 

One final remark before concluding is that detailed anal- 
yses of field statistics would benefit from estimates for (Tpsr) 
in (9), corresponding to the value J^jj in (10). This would 
allow an absolute scale for E' and /' to be estimated accu- 
rately and aid in interpreting detailed fits at low E' of the 
field statistics (Paper III). 



10 CONCLUSIONS 

Analysis of rapidly-sampled, coherently dedispersed data 
show for the first time that Vela pulsar's variability cor- 
responds to field statistics that are (i) well defined and 
(ii) evolve smoothly from Gaussian intensity statistics off- 
pulse to power-law and lognormal statistics on-pulse. Since 
theories for wave growth in inhomogeneous plasmas pre- 
dict the field statistics, these observations allow the source 
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plasma and emission physics to be probed. Detailed single- 
component fits to the observed field statistics confirm that 
the off-pulse variability corresponds to Gaussian intensity 
statistics, consistent with superposition of multiple incoher- 
ent signals and/or scattering, while the lognormal statis- 
tics observed near the peak of Vela pulsar's average profile 
are consistent with the predictions of SGT for a purely lin- 
ear system near marginal stability. The simplest interpre- 
tations are that the Vela pulsar's on-pulse variability is a 
direct manifestation of an SGT state and that only linear 
emission mechanisms (either direct or indirect) are viable. 
This argues against some nonlinear theories for pulsar radio 
emissions. At a deeper level nonlinear mechanisms arc not 
ruled out but arc strongly constrained: viable mechanisms 
must produce lognormal statistics when suitably ensemble- 
averaged. Accordingly mechanisms involving wave collapse 
do not appear viable due to them yielding power-law statis- 
tics. The data suggest that scattering minimally affects the 
field statistics at on-pulsc phases for Vela. Two-component 
fits for Vela's power-law and lognormal domains, described 
elsewhere (Paper III), extend and strengthen these results, 
finding that the observed variability is consistent with log- 
normal statistics whenever the average pulse profile is above 
background. 

These analyses of pulsar emissions represent the first 
applications of SGT to both non-solar-system sources and 
to propagating, free-space radio emissions. This generalizes 
the theory's applicability significantly beyond the relatively 
localized, solar system, plasma waves considered previously. 
While SGT thus applies in all 8 analyses performed by us to 
date, other field statistics are sometimes observed for other 
wave phenomena. Analysis of field statistics is thus a power- 
ful tool for understanding source variability and constraining 
emission mechanisms and source characteristics. This new 
window into the source physics is likely to be widely useful 
for coherent astrophysical and solar system radio emissions, 
as already found for plasma waves in space. 
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